Tracing the Evolution of Temperature in Near Fermi Energy Heavy Ion Collisions 
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The kinetic energy variation of emitted light clusters has been employed as a clock to explore the 
time evolution of the temperature for thermalizing composite systems produced in the reactions of 
26A, 35A and 47A MeV "^"^Zn with '^^^Ni, ^^Mo and ^^'^Au. For each system investigated, the double 
isotope ratio temperature curve exhibits a high maximum apparent temperature, in the range of 10- 
25 MeV, at high ejectile velocity. These maximum values increase with increasing projectile energy 
and decrease with increasing target mass. The time at which the maximum in the temperature 
curve is reached ranges from 80 to 130 fm/c after contact. For each different target, the subsequent 
cooling curves for all three projectile energies are quite similar. Temperatures comparable to those 
of limiting temperature systematics are reached 30 to 40 fm/c after the times corresponding to 
the maxima, at a time when AMD-V transport model calculations predict entry into the final 
evaporative or fragmentation stage of de-excitation of the hot composite systems. Evidence for the 
establishment of thermal and chemical equilibrium is discussed. 

PACS numbers: 25.70.Pq, 24.60.Ky, 05.70.Jk 

Keywords: Liquid gas phase transition, critical fluctuation, fragment topological structure 



INTRODUCTION 

The light particle emission which occurs during violent 
collisions of two heavy nuclei carries essential information 
on the early dynamics and on the degree of equilibration 
at each stage of the reaction. To obtain more specific 
information on the reaction dynamics and on the ther- 
mal evolution of multi-fragmenting systems produced in 
near Fermi energy collisions P, yj^ y] , we have recently 
focused on efforts to investigate the nucleon and light 
cluster emission which occurs prior to disassembly as the 
system thermalizes and equilibrates. In some previous 
works, we have employed coalescence model analyses to 
probe the early dynamic evolution of the reacting sys- 
tem 1^1^. In this paper, we report on the use of similar 
techniques to explore the temperature evolution of hot 
nuclei produced in a series of reactions of 26A, 35A and 
47A MeV 64Zn projectiles with ^SNi, '^^Uo and ^^^Au 
target nuclei using a combined 47r Charged Particle - 47r 



Neutron Ball detection system. The data provide ex- 
perimental evidence for an initial rapid thermalization 
of the incident energy, into a participant matter subsys- 
tem. The double isotope ratio temperature first rises to 
a maximum, then decreases as further particle emission, 
expansion and diffusion of the excitation energy into the 
remainder of the composite system occurs. A close corre- 
lation between the peak temperatures and spectral slope 
temperatures for early emitted particles is also observed, 
suggesting local chemical and thermal equilibration for 
this early emitting system. Temperatures which are com- 
parable to those of limiting temperature systematics 
are reached about 30 to 40 fm/c after the peak temper- 
atures are observed. 
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EXPERIMENTAL DETAILS 

The reactions of 26A, 35A and 47A MeV ^''Zn projec- 
tiles with ^^Ni, ^^Mo and -'^^^Au target nuclei were stud- 
ied at the K-500 Super-Conducting Cyclotron at Texas 
A&M University, using the 47r detector array NIMROD. 
NIMROD consists of a 166 segment charged particle ar- 
ray set inside a neutron ball. The charged particle array 
is arranged in 12 concentric rings around the beam axis. 
The eight forward rings have the same geometrical design 
as the INDRA detector, but have less granularity In 
those rings, the individual segments are fronted by ioniza- 
tion chambers (IC) filled with 30 Torr of CF4 gas. Front 
and back windows were made of 2.0 /im aluminized Mylar 
foil. In each of these forward rings, two of the segments 
have two Si detectors( 150 and 500 /im thick) between the 
IC and Csl detectors (super telescopes) and three have 
one Si detector (300 ^m thick). Each super telescope is 
further divided into two sections. The Csl detectors are 
10 cm thick Tl doped crystals read by photomultiplier 
tubes. For these detectors, a pulse shape discrimination 
method is employed to identify light particles 0. In all 
telescopes particles are, identified in atomic number. In 
the super telescopes, all isotopes with atomic number Z 
< 10 are clearly identified. 

The energy calibration of the Si detectors was carried 
out using both a ^^^Th alpha particle source and the ob- 
served punch through energies of identified reaction prod- 
ucts. The punch through energies were calculated using 
a range-energy table 9]. Since the energy losses of the 
lighter particles, in particular the high energy Hydrogen 
isotopes, are rather small in the Si detectors, evaluation 
of the energy deposited in the Csl crystal from the energy 
loss in the Si detectors requires special care for higher 
energy particles. Therefore, an additional energy cali- 
bration was performed to measure energy spectra from 
the reaction ^''Zn + ^^Mo at 47A MeV. In this run. Si 
detectors of thicknesses 1mm, backed by Csl detectors of 
three different lengths (1cm, 3cm and 5cm), were used 
to measure the inclusive energy spectra of light charged 
particles. The energy spectra were measured at all an- 
gles corresponding to those of the 12 rings of NIMROD. 
The combination of thicker Si AE detectors and observa- 
tion of high energy punch-through points for the particles 
which traversed these thinner Csl detectors allowed us to 
determine the energy spectra with a high degree of con- 
fidence. We then used the '^'^Zn + ^^Mo at 47A MeV as 
a standard reaction to determine the Csl energy calibra- 
tions for all other runs. 

Neutron multiplicity was measured with the An neu- 
tron detector surrounding the charged particle array. 
This detector, a neutron calorimeter filled with Gadolin- 
ium doped pseudocumene, consists of two hemispherical 
end caps and a cylindrical mid-section. The mid-section 
is divided into four separate 90 degree quadrants. The 



hemispheres are 150cm in diameter with beam pipe holes 
in the center and they are upstream and downstream of 
the charged particle array. Thermalization and capture 
of emitted neutrons in the ball leads to scintillation which 
is observed with phototubes providing event by event de- 
terminations of neutron multiplicity but little informa- 
tion on neutron energies and angular distributions. Fur- 
ther details on the detection system, energy calibrations 
and neutron ball efficiency may be found in reference |0| . 

During the experiment, data were taken employing two 
different trigger modes. One was a minimum bias trig- 
ger in which at least one of the Csl detectors detected a 
particle. The other was a high multiplicity trigger which 
required detected particles in 3-5 Csl detectors (depend- 
ing upon the reaction studied). 



DATA ANALYSIS 

An inspection of the two dimensional arrays depict- 
ing the detected correlation between charged particle 
multiplicity and neutron multiplicity in NIMROD (not 
shown) , reveals a distinct correlation in which increasing 
charged particle multiplicity is associated with increas- 
ing neutron multiplicity. Although there are significant 
fluctuations reflecting both the competition between dif- 
ferent decay modes and the neutron detection efficiencies, 
these correlations provide a means for selecting the more 
violent collisions. For the analysis reported in this paper, 
we have selected events corresponding to the largest ob- 
served neutron and charged particle multiplicities. This 
selection corresponds to the 10% of the minimum bias 
trigger events with the highest total multiplicity and em- 
phasizes the lower impact parameter collisions. We refer 
to these events as violent collisions. Many of the tech- 
niques applied in this analysis have been discussed pre- 
viously in greater detail in references 0, 0, Q| . Only a 
brief summary of these is included in the present work. 



Moving source analysis 

A common technique to characterize light particle 
emission in this energy range is to fit the observed spectra 
assuming contributions from three sources; a projectile- 
like (PLF)source, an intermediate velocity (IV)source 
and a target-like (TLF)source. For asymmetric collisions, 
such fits typically exhibit a PLF source dominance local- 
ized at high rapidity, an IV source dominance at mid- 
rapidity and TLF source emission localized at low rapid- 
ity 0, IM ll^ I13I . In the present work, except for the 
most forward detector rings, the data are dominated by 
particles associated with the IV and TLF sources and a 
good reproduction of the observed spectra is achieved. In 
this analysis, the source velocities, temperatures, particle 
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FIG. 1: Slope temperatures from 3-source fits to the experimen- 
tal spectra. Open circles represent apparent temperatures for the 
emission from the IV source. Filled circles represent apparent tem- 
peratures for the emission from the TLF source. 



multiplicities and emission barriers for the three different 
sources were the parameters searched. 

In Figure 1, the slope temperature parameters for emis- 
sion from the IV and TLF sources, derived from the 
fits, are shown for the nine reactions studied. For both 
sources the measured spectra result from a summation 
of the spectra of particles emitted over a range of time. 
Thus the observed slope temperature values are affected 
by the relative emission probabilities over that time pe- 
riod. The IV source slope temperatures for p, d, t, '^He 
and "^He range from T ~ 7 to 17 MeV for the 

different systems studied. The temperatures from the 
particles with A < 3 are quite similar. They follow the 
trends of earlier reported values for pre-equilibrium emis- 
sion at such projectile energies The IV source 
temperatures derived from the alpha spectra are typi- 
cally lower than those measured for the other particles. 
These softer slopes for the alpha particles appear to re- 
flect larger relative contribution of lower energy alphas 
which are attributed to the IV source in the fitting pro- 
cedure. The slope parameters for the TLF sources are 
much lower, in the range of T~ 2 — 6 MeV. For this 
source, the apparent temperatures for alpha emission are 
the highest. Such an effect has previously been noted and 
attributed to the relatively higher emission probabilities 
for alpha particles in the early stage of the evaporation 
cascade [la, 113 ■ 



Temperature Determinations 

Given the continuous dynamic evolution of the system, 
source fit parameters should be considered as providing 
only a schematic picture of the emission process. Never- 



theless, the information derived can be very instructive. 
We have employed such analyses to estimate the mul- 
tiplicities and energy removed at various stages of the 
reaction. To follow the time evolution of the system in 
more detail, a more sophisticated analysis of the par ticle 
emission is necessary. Both theoretical models [Til 
and experiments |2oL l2l| indicate that the early colli- 
sion dynamics leads to a correlation between emission 
time and energy for the early emitted particles. This 
correlation can be exploited to follow the time evolution 
of the system. We have previously used this correlation 
in coalescence model studies of several systems 0, 0, ^3 ■ 
For the present systems, extensive calculations have been 
made employing the AMD-V model of Ono Many 
of the results of the AMD-V calculations for the present 
systems are compared to experimental observables in ref- 
erence [lOj. In the present work, we have employed the 
calculated correlation between particle energy and time 
predicted by those AMD-V calculations to calibrate the 
emission time scales for these reaction systems. 

At intermediate energies, the observed spectral slope 
parameters derived from the source fits are not adequate 
as measures of the temperature evolution, as the observed 
spectra are convolutions of the spectra at different emis- 
sion times and excitation energies and include high en- 
ergy particles which are emitted prior to the achievement 
of thermal equilibrium. For a system at chemical and 
thermal equilibrium at a suitably low density, Albergo 
et al. [23 | have shown that the temperature of the emit- 
ting system can be derived directly from the first chance 
emission 

double isotope yield ratios of two adjacent isotopes of 
two different elements. In a more recent work by Kolomi- 
ets et al. [l^l , essentially the same result is derived when 
only thermal equilibrium is initially assumed. Therefore, 
to characterize the temperature at a particular emission 
time, we have employed double isotope yield ratio mea- 
surements. In the case of strong system evolution, double 
isotope yield ratio temperatures derived from integrated 
yields are certainly suspect if the isotopes being utilized 
are in fact produced at very different times or by dif- 
ferent mechanisms. However, if chemical equilibrium is 
achieved and the particles corresponding to particular 
emission times can be selected, derivations of double iso- 
tope yield ratio temperatures as a function of emission 
time should allow us to follow the temperature evolution 
of the system. 

To focus on the early evolution of the temperature, we 
have first selected clusters observed at mid-rapidity, i.e., 
specifically those detected at angles between 70° and 80° 
in the IV source frame. In this way, we attempt to isolate 
the emission associated with the IV source which occurs 
during the thermalization stage of the reaction 0, 0, 0| . 
We have then made double isotope ratio temperature de- 
terminations as a function of ejectile velocity in the IV 
frame. The velocities used are the "surface velocities" 
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of the emitted particles. The surface velocity, Vsurf, is 
defined as the velocity of an emitted species at the nu- 
clear surface, prior to acceleration in the Coulomb field 
ilia. Vsurf J is obtained in our analysis by subtraction 
of the Coulomb barrier energy derived from the source 
fits. Since the early emitted light particle energies are 
strongly correlated with emission times, and evaporative 
or secondary emission contribute to the spectra primarily 
at the lower kinetic energies, the yields of higher energy 
particles should be relatively uncontaminated by later 
emission processes. This is an important advantage in 
such double isotope yield ratio determinations. At mid- 
rapidity there is little contribution from the PLF source. 
There is, however some observed contribution from the 
TLF source at low velocities in the IV source frame. Since 
the three source fits are only approximations to the emis- 
sion from the continuously evolving system, particles in 
this low velocity range may be viewed alternatively as 
the last particles emitted from the IV source or the earli- 
est from the TLF source. In the following analysis yields 
assigned to the TLF source have been subtracted from 
the experimental yields. 

The temperatures employed are Thhc, derived from 
the yields of d, t, ^He and ^He clusters. For particles 
emitted from a single source of temperature, T, and hav- 
ing a volume Maxwellian spectrum (y^exp — e/T), where 
e is the particle energy, the HHc double isotope yield ra- 
tio evaluated for particles of equal Vsurf ? is y i times the 



surf 7 

ratio derived from either the integrated particle yields or 
the yields at a given energy above the barrier 22] . Thus 



T = 



14.3 



(\/f(l-59i?..„.,; 



In 



(1) 



where the constants 14.3 and 1.59 reflect binding energy, 
spin, masses and mass differences of the ejectiles. If Y 
represents a cluster yield, R{Vsurf) = YdXiHe/YtYiH^ for 
clusters with the same surface velocity. 



Calibration of Timescales 

To calibrate the time-scale associated with our data. we 
have employed results of the AMD- V calculations [lOj . 
In the AMD-V calculations, the particle emission starts 
at times near 50 fm/c, and varies depending on projectile 
velocity and entrance channel masses. From that point, 
the calculations predict an initial rapid decrease of the 
average kinetic energies of the emitted particles with in- 
creasing time of emission. This is followed by a much 
slower rate of decrease at later times, 140 to 200 fni/c 
for the systems studied. The time of the transition from 
rapid to relatively slow kinetic energy depends on projec- 
tile velocity and entrance channel mass. Such trends are 
typically observed in transport model calculations. For 
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FIG. 2: Correlation of average emission time with surface velocity 
for early emitted nucleons as calculated by the AMD-V code. Solid 
symbols depict the results for the nine different reactions as labeled. 



the four systems studied, we have derived from AMD- 
V calculations, which covered the range of time from the 
time of contact up to to 300 fm/c, the correlation between 
average emission times of emitted neutrons and protons 
and their energies. For each reaction the resultant time- 
surface velocity relationship is represented in Figure 2. 
As is seen there, the calculations indicate a near linear 
decrease of Y surf , from near projectile velocity into the 
3 to 3.5 cm/ns range as the average emission times in- 
crease from ~ 50 fm/c to 150 fm/c. In the following we 
have employed these correlations to determine the aver- 
age emission times corresponding to particular observed 
values oiYsurf in the experiment. However, since below 
3.5 cm/ns the experimental data contain large contribu- 
tions from TLF evaporative emission, the sensitivity of 
the emission energy to time is significantly reduced, and 
we do not attempt to assign emission times for particles 
with surface velocities below 3.5 cm/ns. 



RESULTS AND DISCUSSION 

We present, in Figure 3, experimental results for the 
double isotope ratio temperatures as a function of ve- 
locity in the IV frame. For each system investigated, 
the double isotope ratio temperature determination ex- 
hibits a high maximum temperature in the range of 10- 
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27 MeV. These maximum temperatures are much higher 
than the hmiting temperatures determined from caloric 
curve measurements in similar reactions . In each case, 
the apparent temperature decreases monotonically on ei- 
ther side of this maximum. The AMD-V model calcula- 
tions 0, 01 indicate a significant slowing in the rate of 
kinetic energy change in the 3-3.5 cm/ns velocity range. 
The shaded vertical bars in Figure 3 indicate that ve- 
locity value. These points signal the end of the IV (or 
pre- equilibrium) emission stages. At lower velocities, 
the slower nuclear de-excitation modes evaporation, fis- 
sion and/or fragmentation determine the properties of 
the ejectile spectra. At these lower velocities, values of 
T^HHe are 3 to 4 MeV, similar to those spectral integrated 
values seen in other experiments [j^ IzS . I25I Efl . |^. 
These values are also very similar to Thhe temperatures 
calculated when the sequential evaporation code GEM- 
INI is used to simulate the de-excitation of the TLF 
source We take this as further evidence that the 

spectra at these lower velocities still contain contribu- 
tions from late stage evaporation. Temperatures derived 
from the yield ratios in this velocity range require cor- 
rections for secondary decay effects. 

Figure 4 presents the derived Thhc temperatures as a 
function of time. While all nine reactions show a qual- 
itatively similar evolution with time, we now see that, 
for each energy, the time at which the maximum in the 
temperature curve is reached increases with increasing 
target mass and decreases with increasing projectile en- 
ergy. The former observation suggests a longer period 
required for establishment of a thermal and/or chemi- 
cal equilibrium as the total system size increases while 
the latter observation suggests a more rapid thermaliza- 
tion of the initial projectile energy for the initially faster 
projectiles. The dynamic transport calculations indicate 
that the condition of thermal equilibrium of the whole 
system is not yet established as the earliest ejectiles are 
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emitted. Non equilibrium effects may be most evident in 
the particularly high apparent peak temperatures for the 
35A and 47A MeV ^'^Zn + ^^Ni case. The AMD-V calcu- 
lations for the different systems predict a higher degree 
of transparency in those reactions 10]. 

Figure 4 also indicates that the time of entry into the 
evaporation or disassembly stage, taken to be the time 
corresponding to surface velocities of 3.5 cm/ns, increases 
with target mass, from ~ 135 fm/c for the Ni target to 
~ 165 fm/c for the Au target. We note that at such 
times the temperatures are very similar to the limiting 
temperatures derived from a systematic investigation of 
caloric curve measurements 6]. Except for the ®^Zn 
+ ^^Ni reactions at 35A and 47A MeV, the time for the 
initial cooling stage, i.e., the time difference between that 
corresponding to the maximum in the temperature and 
that corresponding to the start of the evaporation stage, 
is 30 to 40 fm/c. 



Interpretation of Temperature Evolution Curves 

The experimental temperature curves show an initial 
rise to a maximum and a subsequent decline. It is tempt- 
ing to interpret the initial rise as reflecting the early rate 
of conversion of projectile kinetic energy into thermal en- 
ergy of the composite system. However, it appears more 
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likely that the double isotope temperature is not accu- 
rately reflecting the temperature of the emitting system 
at earlier times. The time required for the establishment 
of chemical equilibrium is presumably longer than that 
required for thermalization. We are not able to separate 
these with the present data. Rather, in this section we 
address the interpretation of the temperatures and the 
degree to which they can be taken as reflecting thermal 
and chemical equilibration, at least locally if not globally. 



Comparison of Thhc with Tsiope from source fits 

If both thermal and chemical equilibrium are achieved 
and the density is not too high |35j, an agreement be- 
tween the thermal temperature and the double isotope 
ratio temperature Thhc can be expected. 

For the TLF source the fit slope parameters, Tgiope pre- 
sented in Figure 1 are normally lower than the latest time 
"^HHe temperatures seen in Figure 4. This is not surpris- 
ing since the source fit can be expected to return only an 
apparent temperature reflecting the entire cooling stage 
of the TLF source. In our previous work 0, 0, Il7| . 
we have found that the slope parameters for the alpha 
particle emission from the TLF source most closely ap- 
proximates the intial thermal temperature of this source, 
reflecting the higher fraction of the alpha emission in the 
earlier part of the de-excitation cascade. The alpha par- 
ticle slope parameters are presented in column 5 of Table 
I. The relationship between these thermal fit apparent 
temperatures and the late time chemical temperatures 
obtained from the TLF isotope ratio temperatures, as 
presented in column 6, is quite reasonable, given the as- 
sociated uncertainties. 

The observed time evolution of Thhc also prompts 
us to inquire about the relationship between the maxi- 
mum temperatures and the slope temperature parame- 
ters which characterize the ejectile kinetic energy spec- 
tra of the IV sources. These latter are also determined 
from the three source fits to the experimental spectra. 
As is seen in Table I and Figure 5, there is actually a 
close agreement between these two values for all but the 
35 and 47A MeV ^^Zn + ^^Ni reactions. The appar- 
ently high values of the peak temperatures for these last 
two reactions were already seen above. As noted there, 
this may reflect non-equilibrium effects resulting from the 
very high degree of transparency for those two reactions. 

At first glance, the good agreement seen for the rest of 
the reactions still seems surprising because the IV spec- 
trum is known to result from a convolution of the ejectile 
emission from the evolving system. Thus, it depends on 
the time-dependent rates of both the ejectile emission 
and the temperature evolution of the system. However, 
if the requirement that the global fit to the IV source 
reproduce the high energy tail of the IV spectrum plays 
a dominant role in determining the overall slope tem- 
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perature parameter of the fit, this result is understand- 
able. In such a case, the agreement seen between this 
thermal slope parameter and the chemical temperature 
determined from the isotope yield ratios for the lowest en- 
ergy 26A MeV ^*Zn + ^^Ni reaction and for the ^''Zn -f 
^^Mo and ^'^Zn + -'^^^Au reactions at all three projectile 
energies provides a strong indication that a simultane- 
ous thermal and chemical equilibrium is achieved in the 
emitting system, at least locally, at times corresponding 
to those of the peak temperatures. 



Model estimates of peali temperature and source 
sizes 

We have attempted to make a first order estimate of 
the maximum temperature to be expected. In the AMD 
calculations for our systems, only a small fraction of the 
available mass and excitation energy has been removed 
from the system by pre-equilibrium emission at the times 
which we associate with the observed maximum temper- 
atures 0. For the reaction being considered, we then 
ignore the mass loss and estimate the maximum thermal 
excitation energy to be E* = Ecm + Q, where Ecm is 
the available Center of Mass energy and Q, the fusion 
reaction Q value. Excitation energies per nucleon of 5.9 
to 11.2 MeV/u are obtained in this way. The results of 
these calculations are presented in column 7 of Table I. At 
the times corresponding to the maximum observed tem- 
perature, the AMD-V calculations show the composite 
systems to have rebounded from an initial small compres- 
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TABLE I: Temperatures and Excitation Energies 



Target 


Projectile 


Tslope{Avg) 


Max Thhb 


T,;ope(TLF) 


THHe{TLF) 


E* 






Energy 


IV 




alpha 


{vsurf = 3.5cm/ns) 




^part—^tot 

P = Po 




(MeV/u) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV/u) 


(MeV) 




26.0 


9.68 ± 0.68 


9.74 ± 0.25 


5.66 ± 0.30 


4.34 ± 0.40 


5.92 


9.49 




35.0 


11.4 ± 1.05 


16.6 ± 0.99 


5.91 ± 1.43 


5.06 ± 0.45 


8.17 


11.1 




47.0 


13.7 ± 1.28 


24.8 ± 2.57 


6.12 ± 0.80 


5.34 ± 0.46 


11.2 


13.0 


^ Mo 


26.0 


9.60 ± 0.97 


10.0 ± 0.42 


4.18 ± 0.50 


4.50 ± 0.40 


5.55 


9.19 




35.0 


11.9 ± 1.44 


13.3 ±0.99 


5.50 ±0.50 


5.74 ±0.54 


7.73 


10.8 




47.0 


15.1 ± 1.73 


17.6 ± 1.13 


5.41 ± 0.30 


5.71 ±0.53 


10.6 


12.7 




26.0 


9.95 ±0.69 


9.02 ±0.18 


5.13 ±0.13 


4.92 ±0.45 


3.96 


7.76 




35.0 


11.8 ±0.81 


11.9 ±0.31 


5.22 ± 0.35 


5.81 ±0.54 


5.62 


9.25 




47.0 


14.1 ±0.78 


15.7 ±0.54 


6.54 ±0.32 


6.53 ±0.60 


7.84 


10.9 



sion and the system density, p, to be at or below normal. 
If p/po = Ij we can estimate the maximum temperature 
which we might expect. Here we assume a uniform nor- 
mal density Fermi gas |29j with a Fermi Energy which 
we take from the interpolation or extrapolation of the 
values reported in reference \?>(\ . Since the system is 
quite excited, we further assume a nucleon effective mass 
of 1 js^. Ignoring the small mass loss expected, we cal- 
culate the temperature, T, from E* = oT^, where a, the 
level density parameter, is determined from the Fermi 
energy. Here a = A/ 15. 2 MeV~^ is used jlOl- The tem- 
peratures thus derived are presented in Table I, column 
8. While the calculated temperatures reported in column 
8 are indeed significantly higher than limiting termpera- 
tures of caloric curve measurements, they are not as high 
as the observed maximum temperatures in column 4. 

The calculated values in column 8 follow from an as- 
sumed thermalization of the entire system. In a re- 
cent investigation, Sood and Puri have employed a QMD 
transport model to calculate the maximum and aver- 
age temperatures and densities achieved in symmetric 
or near symmetric heavy ion collisions at ^bai, the bal- 
ance energy corresponding to the transition from posi- 
tive to negative flow [2^ . Their calculation of the max- 
imum temperature, based upon a local density approx- 
imation for the matter contained in a sphere of 2 fm 
radius around the center of mass of the system, clearly 
indicates that the entire system is not equilibrated at 
the early times. In general, the available center of mass 
energies in the calculations of reference |2^ are some- 
what higher than those of our reactions with 47A MeV 
^^Zn and the calculations are made for varying impact 
parameters. Nevertheless, the maximum double isotope 
ratio temperatures derived in the present experimental 
study are quite comparable to those reported in reference 
lU This is particularly noteworthy because the calcula- 
tions by Sood and Puri strongly suggest the presence of 
an initial hot, locally equilibrated, participant zone sur- 



rounded by colder spectator matter and the derived tem- 
peratures are interpreted as representing a local thermal 
equilibrium. A similar picture is obtained in the AMD-V 
calculations of reference (T^. Further, such a result is 
consistent with results of earlier experimental studies of 
pre-equilibrium emission which found that the IV spectra 
could be equally well modeled either within the frame- 
work of nucleon-nucleon collision dynamics or as emission 
from a hot thermalized participant zone [l^ lisj . 

If, in fact, the early system consists of both partici- 
pant (nascent fireball |32|) and spectator matter then, 
initially, the available energy may be distributed only 
over a subset of the nucleons. Further, the density for 
this subsystem need not be po- Assuming still that such 
a hot participant zone may be modeled as a uniform den- 
sity Fernii gas allows us to write the more general expres- 
sion mEi 

T ^ \l {Kaip/po)HEx/ Apart), 

where Kq is the inverse level density parameter at nor- 
mal density (in an excited nucleus) and Apart is the mass 
number of the participant zone. Taking Kq to be 15.2 
MeV and the observed peak temperatures from column 
4 of table I allows us to calculate Apart as a function 
of p/po for each system studied at each projectile en- 
ergy. The results are presented in Figure 6 (a)-(c), where 
the calculated Apart, normalized to the projectile mass, 
is plotted against p/po- Some initial compression and 
subsequent expansion is predicted by the AMD-V cal- 
culation. The density at the time of the experimental 
peaking of the temperature is expected to be less than 
normal density. The results of this calculation indicate 
that the number of participant nucleons at the time of 
the peaking of the temperature decreases with increas- 
ing projectile energy and increases with increasing target 
mass. The values of the ratios which would result from 
participation of the total entrance channel masses are in- 
dicated by the horizontal lines in the figure. Most experi- 
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FIG. 6: Calculated values for the correlation of 
p/po consistent with measured peak temperatures. 
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mental values are well below those. These results suggest 
the existence of a hot zone at early times. For the case 
of 47A 64Zn + i^^Au^sults of the AMD-V calculations 
reported in reference[lO| indicate the existence of such a 
zone at early time, with a mass number about twice that 
of the projectile. As seen in Figure 6, the Fermi gas es- 
timate for this system at below normal densities is quite 
close to that. 

While the possibility of emission from a thermalized 
hot zone is one of the possible interpretations which 
has previously been suggested in earlier comparisons of 
dynamic and thermal pictures of pre-equilibrium emis- 
sion in similar collisions 0, 0|, in the present study 
it is inferred not from the slope parameters of the pre- 
equilibrium source but rather from the peak value of the 
temperature, Thhc- The double isotope ratio temper- 
ature measurements assume chemical equilibrium. We 
return to this point in the following section. 



SUMMARY AND CONCLUSIONS 



curve exhibits a high maximum apparent temperature, in 
the range of 10-25 MeV. The maximum values increase 
with increasing projectile energy and decrease with in- 
creasing target mass. These maxima occur at times from 
80 to 130 fm/c after the nuclei contact. They are much 
higher than the limiting temperatures determined from 
caloric curve measurements in similar reactions 0. For 
most of the reactions studied, a close correlation is ob- 
served between the peak temperatures for early emit- 
ted particles, obtained from double isotope ratios, and 
the spectral slope temperatures for the pre-equilibrium 
(IV) source. The data indicate that at least a local ther- 
mal and chemical equilibrium is established during these 
times. After peaking, the temperatures decrease rapidly, 
apparently reflecting particle emission, diffusion of the 
excitation energy into the remaining system and expan- 
sion. For each individual target nucleus, the later por- 
tions of the cooling curves for all three projectile energies 
are very similar, indicating that hot nuclei with similar 
properties are produced. Temperatures comparable to 
those derived from limiting temperature systematics are 
reached 30 to 40 fm/c after the times corresponding to 
the maxima, at the times when AMD-V transport model 
calculations predict entry into the final evaporative or 
fragmentation decay of the hot composite system. 

As a final comment, we note that the present data 
suggest that if the Z=l and Z=2 ejected light particles 
are taken to represent the gaseous phase, as is usually 
assumed, the reaction dynamics of central collisions may 
itself lead to a natural situation in which the gas is not 
outside the liquid matter but initially confined inside the 
liquid matter, perhaps facilitating the establishment of 
a liquid-gas equilibrium. This could somewhat mitigate 
arguments that the concept of establishment of a liquid 
gas equilibrium is not tenable in a nuclear collision as 
there is no container to constrain the gas. Of course, 
more recently it has also been argued, on the basis of 
detailed balance, that the actual physical equilibrium is 
not necessary |2a|- 
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The kinetic energy variation of emitted light clusters 
has been employed as a clock to explore the time depen- 
dence of the temperature evolution of thermalizing com- 
posite systems produced in the reactions of 26A, 35A and 
47A MeV ^^Zn with ^^Ni, ^^Mo and ^^^Au. For each sys- 
tem investigated, the double isotope ratio temperature 
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